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THE  CALCULATION  OF  THE  FLOW  FIELD 
ABOUT  A  BLUNTED  9°  CONE 

SECTION  I 

CHARACTERISTICS  CALCULATIONS 


By  E.  Kennedy,  A.  Fields,  andM.Seidman 


SUMMARY 


The  inviscid  fLow  fields  about  a  hemispherically  blunted  cone 
(9°  half-angle)  have  been  calculated  at  five  hypersonic  Mach  numbers. 
The  calculation  is  begun  in  the  subsonic  region  by  the  method  of 
Vaglio-Laurin  and  Ferri  and  extended  into  the  near  wake  by  the 
method  of  characteristics.  The  air  was  taken  to  be  a  mixture  of 
perfect  gases,  either  chemically  frozen  or  in  thermodynamic 
equilibrium.  Results  from  these  calculations  were  used  as  initial 
data  for  rate  and  wake  calculations  presented  in  Sections  II  and  III 


of  this  report. 
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I.  DISCUSSION 

In  order  to  determine  flight  conditions  for  the  flow  field  calculations,  a 

reentry  trajectory  beginning  at  22,000  fps,  300,000  feet  altitude  and  flight 
o 

path  angle  of  79  was  calculated  on  the  Bendix  G-15  Computer,  With  these 
initial  conditions  and  with  W/C^A  constant  at  30  Ib^ft2  ,  a  ballistic  trajectory 
passes  through  the  following  points  in  the  h-v  plane. 


h 

V 

250,000 

22,550  fps 

200,000 

22,475  fps 

150,000 

21, 600  fps 

100,000 

15,875  fps 

80,000 

9,150  fps 

Some  work  of  Bloom  and  Steiger  (summarized  in  Reference  1)  indicated  that 
at  the  three  higher  altitudes  the  flow  would  be  chemically  frozen  downstream  of  i 
the  nose  region  while  at  the  two  lower  altitudes  it  would  be  closer  to  equilibrium. 
Accordingly,  at  the  lower  altitudes  both  equilibrium  and  frozen  calculations  were 
run  and  onLy  the  frozen  case  at  150,  000  feet  and  above.  The  details  of  the  cal¬ 
culation  are  given  by  Vaglio-Laurin  and  Bloom  in  Reference  1  and  need  not  be 
repeated  here.  Some  difficulty  was  experienced  in  beginning  the  frozen  calculation 
from  the  equilibrium  calculation  in  the  subsonic  and  transonic  region.  This 
difficulty  was  traced  to  the  difference  between  the  frozen  and  equilibrium  speeds  of 
sound  and  resulted  in  subsonic  velocities  in  the  frozen  flow  calculations.  This,  of 
course,  vioLates  a  fundamental  requirement  for  the  existence  of  real  characteristics 
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and  stops  the  calculation.  The  solution  was  to  begin  the  characteristics 
calculation  a  little  further  from  the  sonic  line  so  that  the  reduced  "frozen" 

Mach  number  was  still  greater  than  one.  This  allowed  the  calculation  to 
proceed  to  the  end  of  the  body  without  interruption.  In  order  to  develop  the 
flow  field  about  the  base  of  the  body  and  in  the  near  wake,  it  was  necessary 
to  estimate  the  base  pressure  and  the  position  of  the  wake  shock.  The  base 
pressure  was  estimated  by  the  method  of  Chapman  (2)  and,  in  general,  was 
several  times  the  ambient  pressure.  The  wake  shock  began  at  a  point  on  the 
constant  pressure  streamline  one-half  the  base  radius  from  the  axis.  These 
admittedly  approximate  techniques  enable  one  to  bridge  the  base  region  where  the 
flow  processes  are  fundamentally  dissipative  and  are  not  ye^  to  well  understood. 

The  calculation  was  then  extended  downstream  for  about  20  more  radii  where  the 
pressure  on  the  axis  (more  correctly,  on  the  constant  diameter  core  region 
postulated  in  the  near  wake)  approached  that  of  the  atmosphere.  It  should  be 

pointed  out  that  the  pressure  near  the  bow  shock  at  this  station  was  still  considerably 
above  the  ambient. 

Sufficient  data  were  now  available  to  trace  individual  streamlines  from 
the  bow  shock  to  the  beginning  of  the  far  wake.  Four  such  streamlines  for  each 
flight  condition  are  shown  in  Figures  1  through  7.  With  the  pressure  distribution 
along  these  streamlines  given  by  the  characteristics  analysis,  the  rate  calculations 
of  Section  II  were  performed  and  the  results  used  to  start  the  wake  calculations  of 


Section  III. 
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Aside  from  the  rate  and  wake  calculations,  the  distribution  of  flow 
properties  about  the  body  itself  is  not  without  interest.  Figures  8  through  11 
show  profiles  of  thermodynamic  data  (  p,  p  ,  and  T)  and  electron  density  at  a 
station  1  1/2  nose  radii  from  the  nose  and  Figures  12  through  15  show  the 
same  data  at  the  base  of  the  cone.  (In  the  case  of  frozen  flow,  the  electron 

density  was  calculated  by  assuming  the  reaction  N  +  O _ NO++  e”  to 

be  the  only  significant  electron-producing  one  and  it  to  be  in  equilibrium  at  the 
state  determined  by  the  frozen  field).  In  Figures  12  through  15  ,  both  the 
inviscid  (i.e.  characteristics)  calculations  and  the  modification  of  these  profiles 
due  to  the  boundary  layer  are  shown.  The  boundary  layer  was  calculated  by  two 
different  methods  according  to  the  altitude.  At  the  three  lower  altitudes,  the 
method  of  Reference  3  with  (3  =  0  was  used.  The  flow  at  the  wall  was  assumed 
to  be  in  equilibrium  at  the  assumed  wall  temperature  (1000°K).  It  has  been 
shown  (Reference  4)  that  for  Lewis  and  Prandtl  numbers  of  one,  the  mass  frac¬ 
tions  and  stagnation  enthalpy  obey  a  Crocco-type  relationship  and  so  the  solutions 
of  Reference  3  are  directly  applicable  once  the  free  stream  or  edge  conditions  are 
known.  The  edge  of  the  boundary  layer  was  determined  by  matching  the  mass  flow 
within  the  boundary  layer  with  that  in  the  inviscid  layer.  (Note  that  this  requires 
an  iterative  procedure  in  order  to  determine  the  edge). 

In  the  case  of  the  higher  altitudes  (200,000  and  250,000  feet)  the  shear 
in  the  inviscid  flow  approaches  that  at  the  wall  in  the  boundary  layer  solution 
(Reference  5),  and  some  modification  of  the  above  procedure  is  necessary. 

Methods  of  calculation  for  this  case  are  available  in  References  6  and  7.  The 
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method  of  Reference  7  with  the  same  watt  conditions  ,i...  equilibrium  flow 
at  1000°K)  as  before  were  used.  In  this  solution  the  velocity  and  shear  are 
matched  in  the  two  solution.  (Viscous  and  nonviscous)  in  order  to  determine  the 
edge  of  the  viscous  layer.  The  Crocco  relationship  was  assumed  to  remain 
valid  in  this  case  and  furnished  the  concentration  and  stagnation  enthalpy  profiles. 

The  triangular  symbols  on  Figures  12  are  the  values  given  by  the 
nonequilibrium  streamline  calculation.  (Since  the  equilibrium  electron  density 
at  80,000  feet  is  rather  low  (~105-  104/cc>.  this  C4S„  missing  from  Fig_ 
ures  8  through  12.) 
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7 


Streamline 

^0, 

**n2 

NO 

250, 

000 

Body 

.000012 

.236 

.  28 

.482 

.0013 

A 

.00002 

.236 

.  194 

.568 

.0019 

B 

.00015 

.235 

.0300 

.  732 

.0037 

C 

.  104 

.118 

.000014 

.  750 

.0275 

D 

.210 

.0215 

-- 

.  759 

.0099 

200, 

000 

Body 

.00003 

.  235 

.  290 

.  472 

.0025 

A 

.00009 

.  235 

.134 

.  626 

.0044 

B 

.0027 

.227 

O 

o 

.  753 

.0131 

C 

.  143 

.  078 

-- 

.  749 

.0302 

D 

.237 

.  763 

_  — 

150, 000 

Body 

.0001 

.233 

.  21 

.  5535 

.0064 

A 

.0004 

.  231 

.0651 

.  694 

.0102 

B 

.0608 

.  153 

.  0002 

.  744 

.0414 

C 

.  1975 

.0173 

.  000012 

.  745 

.0407 

D 

.237 

-- 

.  763 

_ 

100, 000 

Body 

.0076 

.213 

.  00  7 

.  742 

.0300 

A 

.0476 

.1581 

.0009 

.  7341 

.0593 

B 

.1451 

.0  620 

.000046 

.  7362 

.0567 

C 

.237 

-- 

.  763 

D 

.  237 

-- 

-- 

.  763 

80, 000 

Body 

.22 

.0035 

1.8  x  10'? 

.  7515 

.  025 

No  Dissociation  . 
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